Experimental evidence for water formation via ozone 
hydrogenation on dust grains at 10 K 
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ABSTRACT 



The formation of water molecules from the reaction between ozone (O3) and 
D-atoms is studied experimentally for the first time. Ozone is deposited on non- 
porous amorphous solid water ice (H 2 0), and D-atoms are then sent onto the 
sample held at 10 K. HDO molecules are detected during the desorption of the 
whole substrate where isotope mixing takes place, indicating that water synthesis 
has occurred. The efficiency of water formation via hydrogenation of ozone is of 
the same order of magnitude of that found for reactions involving O atoms or 
O2 molecules and exhibits no apparent activation barrier. These experiments 
validate the assumption made by models using ozone as one of the precursors of 
water formation via solid-state chemistry on interstellar dust grains. 

Subject headings: astrochemistry — dust, extinction — ISM: molecules — molecular 
processes — methods: laboratory 
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Introduction 



Solid water is a very abundant material. Indeed, it is believed to be the most abundant 
condensed material in the Universe, thanks to its propensity for remaining as a deposit 
on interstellar dust particles in dense clouds. Amorphous solid water ice, together with 
a wealth of other condensed species, is widely observed spectroscopically via infra-red 
absorption lines present in the spectrum of a field star o r of the embedded object in star 



forming regions (IGibb et al. 



2000 



Boogert et al. 



2004 ) . Icy mantles covering dust grains 



are reputed to form after direct freeze out of gas phase species and after surface reactions 
of atoms and radicals. 

Although water ice is by far the most abundant species in icy mantles, the chemical 
origin of water molecules in interstellar icy mantles is still uncertain. Water molecule 
formation in the gas phase is not efficient e nough to reproduce the observed abundanc es 



in dark clouds, especially in its solid form (IParise et al. 



2005 



Ceccarelli et al. 



20071). 



Therefore, water ice is likely to form directly on the cold interstellar grains and not as a 
condensate after formation in the gas phase. 



A reaction scheme inc 
ago by 



Tielens fc Hagen 



uding water formation on grain surfaces was proposed some years 



(119821 ). They suggested that H 2 formation would be initiated by 
H-atoms reacting with O, O2 and O3, although the O3 + H pathway was considered to be 
the most effective and O2 would play a mere catalytic role. 

Nonetheless, O2 and O3 molecules have not so far been detected in interstellar ices 
on dust grains, although their presence is presumed. Ozone non-detection may be partly 
explained because its strongest u 3 (O-O) vibrational transition band at 1040 cm" 1 (9.6 
jttm) is swamped by the strong infrared absorption band of silicates at 10 /^m which has 
been observed by ISO (Infrared Space Observatory) in the 2.4-25 fim r egion (e.g., spectrum 



toward the protostar W33A, resolution ~ 500-1000, 



Gibb et al. 



2000) 



-4- 



Astrochemical models of dark clouds predict that condensed oxygen is likely to be a 



major component 
O + H pathwa y ( 



Ioppolo et al. 



of apolar dust grai ns within interstellar clouds at 10 K. Recent 



Du 



2008 



ieu et al. 



Matar et al. 



2009) and the O2 + H pathway (IMiyauchi et al. 



y, the 



2008 : 



20081 ) have been explored experimentally. All these 
experiments show that these reactions lead to water formation, even though the efficiencies 
of intermediate reactions are still under debate. Besides O and O2, ozone too can be present 
in icy mantles chemically pr ocessed by irradiation from Galactic cosmic-ray particles and 



internal secondary UV field (ITielens fc Hagen 



Strazzulla fc Johnson 



1982 



Mathis et al. 



1983 



Greenberg 



1984 : 



19911 ). Several laboratory works confirm the efficient production of 



O.s following energetic electron irradiation of pure oxygen ices (IBennett fc Kaiser 



2005 



Sivaraman et al. 



20071 ). and UV, pro ton and ion irradiation o 



realistic analogues o f interstellar ices (jEhrenfreund et al. 



Cooper et al. 



1992 



ice mixtures simu 



Strazzulla et al. 



atmg 



1997: 



20081 ). In solid phase, O2 and CO2 are easily dissociated by UV radiations 



(A < 240 nm) and cosmic rays (c.r.) producing ato mic O in the electronic ground state 
( 3 P) and in the first electronically excited ( X D) state (ISivaraman et al. 1120071 ): 



„ hv.c.r. „ „ 

2 — — ► O + O 



X 



C0 2 CO + o 



(2) 



In the interstellar grain mantles, the 0( 3 P) atoms in the ground state, which are in thermal 
equilibrium with the ice surface, react with O2 molecules to form O3 via the following 
reaction pathway: 

O + 2 — 3 (3) 
The formation of ozone via equation (3) is exothermic (AfH° = —1.48 eV). 
While the main ozone destruction reactions are: 



O3 + H 



OH + O2 



(4) 
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„ hu,c.r. „ „ 

o 3 — > o 2 + o 



(5) 



In addition, in dark clouds hydro gen is mainly in its mole cular form so H-atoms are a rather 



rare reactant with H/H 2 ~ 10 3 (ILi &: Goldsmith 



20031 ). The number density of H-atoms 



is mostly governed by the destruction of H 2 due to cosmic rays. This value, whatever the 
density of the cloud, is about 1 cm -3 . The 0/H 2 ratio remains approximately constant 
(10~ 4 ), thus the number of atomic O, unlike H, i s proportional to the density of the cloud. 



See, for example, Table 1 of 



Caselli et al. 



(120021 ): for a cloud density of 10 4 cm" 3 the H/O 



ratio is ~ 1/0.75 while for a denser cloud with a density of 10 5 cm -3 the H/O ratio is ~ 1/7. 
Therefore, for very dense clouds, O is the most abundant species in atomic form and it can 
accrete on grains and subsequently form 2 and 3 . 



Cuppen fc Herbst 



(120071 ) showed that 



In their model, using a Monte Carlo approach, 
under dense cloud conditions the (O3 + H) pathway is the most efficient route leading 
to the formation of H 2 0. Finally, it should be noted that, in the diverse models, som e 



barriers can be rais e d in order to increase one of the competitive reactions (ILee et al. 



Caselli et al. 



2002 



Parisc 



1978 



20041 ). Yet, such barriers can only be experimentally probed, 



because extrapolations from the gaseous data are quite uncertain. 

This Letter is the first report of water formation via ozone hydrogenation in the series 
of recent experiments that aim to investigate water synthesis in the ISM. We present the 
evidence for the formation of water molecules from the reaction of O3 with D-atoms on 
non-porous amorphous solid water ice under conditions relevant to the dense interstellar 
clouds. 
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2. Experimental procedures 

Our experiments were performed using the FORMOLISM set-up developed to study 
the reaction and the interaction of atoms and molecules on surfaces simulating dust grains 
under interstellar conditions. The apparatus is composed of an ultrahigh-vacuum chamber 
with a base operating pressure of around 10~ 10 mbar. The sample holder, a copper cylinder 
block, is attached to the cold finger of a closed-cycle He cryostat. It can be cooled to 
8 K. The temperature of the sample is measured using a calibrated silicon diode and a 
thermocouple (Cromel-AuFe) clamped on the back of the sample holder. The copper surface 
is covered with 100 layers of non-porous amorphous solid water (np-ASW) by condensation 
of water vapor on the cold sample maintained at 120 K. Water is sprayed through a 



microchanne 



array doser located 2 cm in front of the surface. This ice is amorphous but 



non porous (IKimmel et al. 



details can be obtained in 



2001). The water ice substrate is then cooled to 10 K. More 



Amiaud et al. 



(|2006|). 



Ozone can be synthesized by two methods. One is an ex-situ method: O3 is prepared 
from gaseous 2 introduced in a glass bottle at a pressure of about 30 mbar. The 2 gas 
is then excited by a radiofrequency electric discharge (~ 2 MHz) through a copper coil 
placed around the 2 container. The blue ozone gas produced by the recombination of 
O2 and O atoms during the discharge is condensed on the bottom wall plunged in liquid 
nitrogen (77 K). Residual O2 is then removed by primary pumping. The glass bottle is then 
mounted onto the main set-up. Ozone is introduced into the vacuum chamber via a triply 
differentially pumped beam line aimed at the cold surface covered with H 2 ice. Due to the 
metallic parts of the inlet system on which O3 in part dissociates, the beam is composed of 
70 % ozone and 30 % 2 . However, 2 is easily removed from the surface when the sample 
is held at 50 K. By varying the exposure time of the ozone beam, it is possible to vary the 
initial coverage of the surface. The O3 flux was previously calibrated using temperature 
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programmed desorption (TPD) by determining the exposure time of a pure O2 beam 
required to saturate the O2 monolayer (1 monolayer = 10 15 molecules cm -2 ) on np-ASW 
ice. In the experiments described below we have deposited a maximum of half a monolayer 
of ozone under the same beam conditions used for calibrating the O2 flux. Owing to the 
non-pure O3 beam, however, ozone exposures had a relative accuracy of ~ 20 % and an 
absolute accuracy of ~ 30 %. Nevertheless, we did not see any significant effect due to the 
initial coverage (between 0.1 and 0.5 monolayers ± 30 %). 

An alternative technique of ozone production was also used, namely by depositing an 
atomic O on top of an O2 layer pre-adsorbed on the surface and then producing ozone 
in-situ. The disadvantage of this technique is that the control of the initial amount of ozone 
have to be checked after each exposure. Whichever technique was used to deposit ozone 
onto the surface, we obtained the same final results. 

After depositing ozone, an atomic D-beam aimed at the surface of the sample irradiates 
the O3-H2O film maintained at 10 K. The D-atoms are produced by the dissociation of D2 
molecules in a quartz tube using a Surfatron microwave discharge at 2.45 GHz. The rate 
of dissociation is ~ 60 %. After the irradiation of ~ 5 x 10 15 D-atoms cm -2 , a TPD is 
performed up to 180 K so that all the water ice desorbs. The TPD heating rate is 10 K/min 
from 10 to 180 K. The desorbed species are detected using a quadruple mass spectrometer 
(QMS). 



3. Results and discussion 

The spectrum (a) in Figure Q] shows the TPD curve of mass 32 obtained after the 
exposure of half a layer of O3 on a np-ASW substrate held at 10 K (ex-situ synthesis). 
We observe two peaks: the desorption peak at low temperature (~ 30 K) represents the 
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desorption of O2 molecules. These molecules have been deposited by the O2 fraction in the 
O3 beam, probably due to ozone destruction on the metallic walls of the gas inlet system. 
The high-temperature peak at ~ 60 K is the desorption of ozone. Mass 48 of O3 is also 
simultaneously detected in smaller amount and is always proportional to the signal at mass 
32. In fact, for the most part fragments are detected, as a result of the O3 cracking 
in the head of the QM S. Also, the desorption of ozone from ASW surfaces is known to 



occur at around 60 K (IChaabouni et al 



2000 



Borget et al. 



2001 



Although this may 

be seen as an indirect measure, mass 32 gives information about O2 (below 50 K) and O3 
(above 50 K) desorption. As shown in the spectrum (a) of Fig. 1, O3 desorption is about 
three times larger than that of O2. This indicates that the ozone fraction deposited with 
the ex-situ method is ~ 70 %, the rest being 2 molecules. In order to deposit only ozone 
molecules on the cold surface, we performed the exposure step with the water ice substrate 
held at 50 K. Thus O2 molecules desorbing at 30 K were eliminated. Then we cooled the 
sample to 10 K before proceeding with the TPD. As shown in the spectrum (b) of Fig. 1, 
the TPD peak observed comes from O3 desorbing at 60 K. By using these procedures, we 
are sure that at 10 K only reactions involving ozone and D-atoms will occur. 

For this work, we performed a set of experiments with various D-atom exposures. In 
the first one, D-atoms were sent onto the ozone-free water ice surface to check that no 
reaction occurred between D and the H 2 substrate to produce deuterated water molecules. 
This result is displayed in the mass 19 spectrum (b) of Fig. 2. There is neither HDO nor 
D 2 desorption peak when D-atoms alone ar e depo sited on the ASW ice sub strate at 10 K 



as was previously noticed by iNagaoka et al. 



(|2005|) and 



Dulieu et al. 



(120091). 



In a second experiment we exposed the np-ASW sample held at 10 K to D-atoms 
before ozone deposition. In this experiment ozone desorbs as if no atoms had been sent 
to the surface. This is because D-atoms form promptly D 2 on non-porous water ice (e.g. 



- 9- 



Congiu et al. 1120091 ) and apparently because ozone does not react with D 2 (see below). 
This proves that hydrogenation reactions take place at around 10 K and not at higher 
temperatures since D-atoms reac t promptly at 10 K and furthermore they will have 



completely desorbed above 13 K (lAmiaud et al. 



20071 ). 



In another experiment, D-atoms were sent onto the np-ASW surface previously exposed 
to 0.5 monolayers of O3. In the spectrum (c) of Figure 1 we firstly observe the disappearance 
of the O3 peak. This suggests that all the ozone molecules have reacted. We also note that 
the reaction D + O3 is very efficient, because the number of D-atoms sent to the surface is 
only 10 times the number of O3 on the surface. Consequently, we can conclude that the 
D + O3 reaction — which is in competition with the D + D reaction in our experiment - 



proceeds on the ice surface with a small energy barri er. An estimation o 



can be made, for example, with the formula used by 
the surface reactions in their model: 

2a 1—- — 

Rreact = V r eXp ( — ^2E a fJ,), 



Cuppen fc Herbst 



the reaction rate 



(120071 ) to simulate 



(6) 



where v T is the attempt frequency for reaction, a the width of the barrier, \x the reduced 
mass and E a . the activation barrier. By ta king E a /k-B = 450 K for the D + O3 reaction 



(ILee et al. 



1978 



Cuppen fc Herbst 



20071 ). we found that this reaction is very rapid 



1 ps) if compared with the diffu sion time of D-atom s on an amorphous water ice 



fa 



iff 



10 ms on porous ASW ice, 



Matar et al. 



20081 ) 



In passing, we would like to point out that, over the last years, we have performed 
experiments on water formation in space starting from O, O2 and now O3. Even though 
we cannot infer absolute reaction rates which are very hard to estimate we can attempt to 
make relative comparisons. It should be noted that the destruction rate of ozone following 
D exposures is as important as the one for O + D and D + O3 reactions within the error bar 
of our experimental data. Hence we conclude that the D + O3 reaction proceeds without a 



-im- 



measurable activation barrier in the solid phase. 

Obtained subsequently to the D + 3 experiment, the spectrum (a) of Figure [2] shows 
a desorption peak of mass 19 at 160 K attributed to the singly deuterated water (HDO). 
This clearly indicates that the reaction of D-atoms with O3 molecules on the ice surface is 
efficient and produces deuterated water molecules detected upon desorption of the whole 
ice film. We also notice that nearly no desorption peak of D2O (mass 20) is detected. This 
i s due to thermally activated H/D exchanges during the phase transition of the water ice 



( iSmith et al. 



19971 ). Therefore, even if D2O is actually formed on the surface, mostly water 



molecules in the form of HDO are likely to desorb and be detected. 

The formation of deuterated water molecules at T s = 10 K from D-atoms and O3 
molecules may then follow the successive reaction pathways: 

3 + D — > OD + 2 (7) 

OD + D — >D 2 (8) 

2 + D^D0 2 (9) 

D0 2 + D^D 2 2 (10) 

D 2 2 + D^D 2 + D (11) 

Whereas at T s greater than 130 K, as mentioned above, the following isotope exchange 
reaction is active: 

H 2 + D 2 — ^2 HDO (12) 

Of course, the reversal of reaction (12) is also possible, but in our case it is statistically 
negligible. 

As for what concerns the intermediate products of the reaction scheme, we performed 
again the D + O3 experiment but with a low exposure of D-atoms (0.01 monolayers) on 
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0.15 monolayers of O3, thus with a dose of D insufficient to go on with the hydrogenation 
reactions leading to water molecules. This experiment shows that O3 molecules are not 
thoroughly consumed and also that O2 molecules appear on the surface indicating that 
reaction (7) proceeds at 10 K and that the formation of D2O can proceed as well from the 
intermediate reaction (9). 

The studied mechanism for water formation from the hydrogenation of ozone might be 
an explanation for the non- detection of O3 in the interstellar ices. In fact, ozone reacts 
rapidly with hydrogen atoms to form water molecules on a water ice surface at 10 K and 
under conditions similar to those found in dense interstellar clouds. Likewise, the lack of a 
clear detection of O2 in interstellar ices may be explained by the high reactivity between O2 
and H-atoms reported in Miyauchi et al. 2008, Ioppolo et al. 2008, and Matar et al. 2008. 

In addition, the OH + H 2 reaction could also be a key surface reaction in dark 
clouds, because Ho (and D2) is present on the surface of the grains in large amounts 



( iKristensen et al. 



2009J). In our experiments too, there is a large amount of molecular 
deuterium because of the incomplete dissociation of D 2 in the D beam. At the moment it 
is not possible to us to probe explicitly the OH + H 2 reaction, but we can investigate the 
role of D 2 , by testing the O3 + D 2 reaction. In order to do that, we have deposited a very 
large amount of D2 (~ 10 17 cm -2 ) on np-ASW ice held at 10 K and previously dosed with 
O3. In this experiment, we did not detect any HDO, even when the surface was held at 
50 K during D 2 deposition in order to overcome a possible activation barrier. Therefore, 
we can fairly conclude that the D 2 + O3 reaction does not proceed under our experimental 
conditions. 

Unlike the O3 + D 2 reaction, O3 + D shows no barrier under our experimental 
conditions and this result shows clearly that ozone is a likely precurso r to water formation 



on a water ice layer at 10 K as proposed by 



Cuppen fc Herbst 



(120071 ). This pathway is to 
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be preferred in environments where the number of O atoms exceeds the number of H atoms, 
which is supposed t o be the case in dark clouds with a density higher than 10 4 cm -3 



(jCaselli et al. 



20021 ). Unfortunately, our present data do not allow us to disentangle the 
different intermediate reactions that lead to the formation of water. In order to investigate 
in more detail the kinetics of water formation at 10 K via ozone, we are developing an 
additional IR diagnostic on FORMOLISM . This will allow us to see the appearance of some 
of the intermediate reactants and to provide accurate quantitative results. 



4. Conclusions 

In this paper we have shown that it is possible to form efficiently water molecules from 
the reaction of O3 with hydrogen atoms on non-porous water ice. This reaction exhibits 
no activation barrier at 10 K, and is about as efficient as the pathways having O or O2 as 
precursors. This experimental result corroborates theoretical studies that consider ozone as 
one of the major actors of water formation in dark clouds. 
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well as the strong financial support from the Conseil Regional d'lle de France (SESAME 
program I-07-597/R), the Agence Nationale de la Recherche (ref. 07-0129) and the Conseil 
General du Val d'Oise. 

Henda.Chaabouni@u-cergy.fr. 



-13- 
REFERENCES 

Amiaud L., Fillion J. H., Baouche S., Dulieu F., Momeni A., Lemaire J. L. 2006, J. Chem. 
Phys, 124, 94702 

Amiaud, L., Dulieu, F., Fillion. J.-H., Momeni, A., Lemaire, J. L. 2007, J. Chem. Phys, 
127, 144709 

Bennett, C. J., k Kaiser, R. I. 2005, ApJ, 635, 1362 

Boogert, A. C. A., Pontoppidan, K. M., Lahuis, F., J0rgensen, J. K., Augereau, J.-C, 
Blake, G. A., Brooke, T. Y., Brown, J., Dullemond, C. P., Evans, II, N. J., Geers, 
V., Hogerheijde, M. R., Kessler-Silacci, J., Knez, C, Morris, P., Noriega- Crespo, A., 
Schoier, F. L., van Dishoeck, E. F., Allen, L. E., Harvey, P. M., Koerner, D. W., 
Mundy, L. G., Myers, P. C, Padgett, D. L., Sargent, A. I. & Stapelfeldt, K. R. 2004, 
ApJS, 154, 359 

Borget, F., Chiavassa, T., Allouche, A., Aycard, J. P. 2001, J. Phys. Chem. B, 105, 449 

Caselli, P., Stantcheva, T., Shalabiea, O. M., Shematovich, V. I., & Herbst, E. 2002, 
Planetary and Space Science, 50, 1257 

Ceccarelli, C, Caselli, P., Herbst, E., Tielens, A. G. G. M., & Caux, E. 2007, in Protostars 
and Planets, ed. V. B. Reipurth, D. Jewitt, & K. Keil (U. Arizona Press), 47 

Chaabouni, H., Schriver-Mazzuoli, L., & Schriver, A. 2000, J. Phys. Chem. A, 104, 6962 

Congiu, E., Matar, E., Kristensen, L., Dulieu, F and Lemaire, J. L. 2009, MNRAS, 397, L96 

Cooper, P. D., Moore, M. H., and Hudson, R. L. 2008, Icarus, 194, 379 

Cuppen, H.M., & Herbst, E. 2007 ApJ, 668, 294 



-14- 

Dartois, E. 2005, Space Science Reviews, 119, 293 

Dulieu, F. Amiaud, L., Congiu, E., Fillion, J. H., Matar, E., Momeni, A., Pirronello V., 
Lemaire J. L., submitted to A&A. 2009:ArXiv 0903.3120vl 

Ehrenfreund, P., Breukers., R. and D'Hendecourt, L. & Greenberg, J. M. 1992, A&A, 260, 
431 

Gibb, E. L., Whittet, D. C. B., Schutte, W. A., Boogert, A. C. A., Chiar, J. E., Ehrenfreund, 
P., Gerakines, P. A., Keane, J. V., Tielens, A. G. G. M., van Dishoeck, E. F. & 
Kerkhof, O. 2000, ApJ, 536, 347 

Greenberg, J. M. 1984, Sci. Am, 250, 124 

Ioppolo, S., Cuppen, H. M., Romanzin, C, van Dishoeck, E. F., & Linnartz, H. 2008, 
Astrophysical Journal, 686, 1474 

Kimmel, G. A., Stevenson, K. P., Dohnalek, Z., Smith, R. S., & Kay, B. D. 2001, J. Chem. 
Phys, 114, 5284 

Kristensen, L. E., Amiaud, L., Dulieu, F., Lemaire, J. L. 2009 submitted to A&A 
Lee, J. H., Michael, J. V., Payne, W. A., & Stief, L. J. 1978, J. Chem. Phys, 69, 350 
Li, D., & Goldsmith, P. F.2003, ApJ, 585, 823 

Matar, E., Congiu, E., Dulieu, F., Momeni, A., & Lemaire, J. L. 2008, A&A Lett, 492, L17 
Mathis, J. S., Mezger, P., & Panagia, N. 1983, A&A, 128, 212 

Miyauchi, N., Hidaka, H., Chigai, T., Nagaoka, A., Watanabe, N., & Kouchi, A. 2008, 
Chemical Physics Letters, 456, 27 

Nagaoka, A., Watanabe, N., & Kouchi, A. 2005, ApJ, 624, L29 



- 15 - 

Parise, B., Ceccarelli, C, & Maret, S. 2005, A&A, 441, 171 

Parise, B. 2004, PhD Thesis, CESR (Toulouse, France) 

Pontopiddan, K. M., van Dishoeck, E.F., Dartois, E. 2004, A&A, 426, 925 

Schmitt, B., de Bergh, C, Festou, M (Eds.). 1998, Solar System Ices, Kluwer, Dordrecht 

Sivaraman, B., Jamieson, C. S., Mason, N. J., & Kaiser, R. I. 2007, ApJ, 669, 1414 

Smith, S. R., Huang, C, & Kay , B. D., 1997, J. Phys. Chem. B 101, 6123 

Strazzulla, G., & Johnson, R. E. 1991, in Comets in the Post-Halley Era, ed. R. L. Newburn, 
Jr., M. Neugebauer, & J. Rahe (Dordrecht: Kluwer), 243 

Strazzulla, G., Brucato, J. R., Palumbo, M. E., & Satorre, M. A. 1997, A&A, 321, 618 

Tielens, A. G. G. M., & Hagen, W. 1982, A&A, 114, 245 



This manuscript was prepared with the AAS IATjrpC macros v5.2. 



-16- 




10 20 30 40 50 60 70 80 90 



Temperature (K) 



Fig. 1. — TPD spectra of mass 32 obtained after depositing ~ 5 x 10 14 molecules cm~ 2 of 
3 on np-ASW: (a) 3 deposited at 10 K, (b) 3 deposited at 50 K and cooled to 10 K, (c) 
3 deposited at 50 K, cooled to 10 K and irradiated with 5 x 10 15 cm -2 of D-atoms. 




Fig. 2. — TPD spectra of mass 19 (HDO) obtained after: (a) deposition of 5 x 10 15 cm~ 2 D- 
atoms on the np-ASW substrate at 10 K previously exposed to 5 x 10 14 molecules cm -2 of 
3 , (b) deposition of 2 x 10 16 cm~ 2 D-atoms on the ozone-free water ice surface at 10 K . 



